States of interacting composite fermions at Landau level fillig v = 2 + 3/8 
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There is increasing experimental evidence for fractional quantum Hall effect at filling factor v — 
2+3/8. Modeling it as a system of composite fermions, we study the problem of interacting composite 
fermions by a number of methods. In our variational study, we consider the Fermi sea, the Pfaffian 
paired state, and bubble and stripe phases of composite fermions, and find that the Fermi sea state 
is favored for a wide range of transverse thickness. However, when we incorporate interactions 
between composite fermions through composite-fermion diagonalization on systems with up to 25 
composite fermions, we find that a gap opens at the Fermi level, suggesting that inter-composite 
fermion interaction can induce fractional quantum Hall effect at i/ = 2 + 3/8. The resulting state is 
seen to be distinct from the PfafRan wave function. 
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I. INTRODUCTION 

The fractional quantum Hall efFeciii (FQHE) is under- 
stood as a consequence of the formation of bound states 
of electrons and quantized vortices, known as composite 
fermions^.. In particular, the FQHE at fractions belong- 

is a manifestation 



2pn±l 



ing to the sequences v = s ± 

of the integral quantum Hall effect (IQHE) of composite 
fermions at the composite-fermion filling ly* = n. In re- 
cent years, the FQHE at fractions not belonging to the 
these sequences has attracted interest, because it can- 
not be explained in terms of a model of noninteracting 
composite fermions, which only exhibit IQHE. In many 
instances, these new fractions can be shown to arise 
from the weak residual interactions between composite 
fermions. For example, the FQHE at z/ = 5/2 is un- 
derstood as a p-wavc Pfafhan-paired state of composite 
fermion^'i, and the FQHE at 4/11 as a fractional QHE of 
composite fermions^'^'^. Recent experiments in very high 
mobility samples have found signatures^ for FQHE at 
z^ = 2 + 3/8, albeit with a tiny gap of a few mK. Al- 
though the evidence is not yet conclusive, the possibility 
of this FQHE is particularly exciting both because it is 
an even denominator fraction and because this fraction 
occurs in the second Landau level, where FQHE is not 
as extensive as in the lowest Landau level. That has mo- 
tivated us to examine various theoretical possibilities at 
this filling factor proceeding with the assumption that the 
second Landau level 3/8 state can be modeled as com- 
posite fermions at filling 3/2. Within this model, we rule 
out several simple variational states, but find an insta- 
bility of the composite fermion Fermi sea into a gapped 
state; while the true nature of this FQHE state is not 
fully understood at present, our studies indicate that it 
is distinct from the usual Pfaffian state. We demonstrate 
the robustness of this state against finite thickness. 

The strongly interacting system of electrons confined 
to a Landau level is described in terms of exotic emer- 
gent particles called composite fermions, which are bound 



states of an electron and an even number (2p) of vor- 
tices of the many-body wave function. The most dra- 
matic consequence of the composite fermion (CF) for- 
mation is that the Berry phase due to the attached 
vortices effectively cancels part of the external mag- 
netic field, producing dynamics governed by a reduced 
field B* = B — 2pp(j)o, where p is the CF density and 
00 = hc/e is the flux quantum. Composite fermions 
form Landau-like levels, called A levels, with their filling 
factor related to the electronic filling factor through the 
expression v — 2~TFTi- ^^^ CF formation correctly ac- 
counts for most of the correlation effects, and a model of 
weakly interacting composite fermions securely explains 
the prominent experimental observations, including the 
FQHE V = ^ 2±i ^^ *^he IQHE of composite fermions^ 
and the compressible liquid at v — l/{2p) as the Fermi 
sea of composite fermionsiS. More subtle structures can 
emerge due to the weak residual interactions between 
composite fermion o^'^'^^i^^'^^i^^'^^i^^ . 

Treating the lowest filled Landau level as inert, the 
problem of our interest is that of interacting electrons 
in the second LL at v^^^ = 3/8. The dimension of the 
Hilbert space here is too large to obtain meaningful re- 
sults from exact diagonalization. We proceed instead by 
modeling v^^' = 3/8 as a state of composite fermions 
at filling factor i^* = 3/2. Assuming that composite 
fermions are fully spin polarized, the state contains a 
fully occupied lowest A level, and a half filled second A 
level. We assume that the lowest A level is inert and work 
with only the composite fermions of the half filled second 
A level in the rest of the paper, denoting their number 
byA^. 

Several states of composite fermions have been con- 
sidered at a half filled Landau level: Fermi sea, paired 
Pfaffian state, stripes, and bubble crystals. When consid- 
ering the first two states at half-filled second A level, the 
^CFs in the second A level capture two additional vortices 
to transform into ^CFs, thereby producing very complex 
mixed structures. (The symbol ^^CF refers to composite 



fermions carrying 2p vortices.) Which, if any, of these 
states is stabihzed in nature depends on the residual in- 
teraction between composite fermions in their second A 
level, which itself is a remnant of the Coulomb interaction 
between electrons occupying the second Landau level. 

It is interesting to note that when composite fermions 
in the second A-level capture an additional pair of vor- 
tices, as is the case for the CF-Fermi sea or the Pfaf- 
fian state, three species of fermions coexist in the system: 
electrons in the lowest LL, ^CFs in the lowest A-level of 
the second LL, and and *CFs in the second A-level; only 
the last will be explicitly considered in our calculations. 
Much or our effort will be toward obtaining the effective 
interaction between them by integrating out the remain- 
ing fermions. 

Our paper is organized as follows. In Sec.|TT]an effective 
inter-CF interaction is derived for composite fermions in 
the second A level of the second Landau level; both pseu- 
dopotential and real-space representations are obtained. 
Using this interaction several variational states are com- 
pared energetically in Sec. IIIII Exact diagonalization re- 
sults in Sec. IIVI confirm the relative advantage of the CF 
Fermi sea state. Then, in Sec. |Vl the residual interaction 
between composite fermions is included perturbatively 
to explore any further instability of the CF Fermi sea 
state. Sec. I VII summarizes the principal conclusions of 
our study. 



II. INTER-COMPOSITE-FERMION 
INTERACTION 

The determination of the inter-CF interaction proceeds 
along several steps. First of all, following standard prac- 
tice, we represent the second LL Coulomb interaction 
(including finite thickness effects) as an effective interac- 
tion in the lowest LL (with zero thickness), for which we 
use the form— 
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(1) 



where the constants B^ , Ci are evaluated by matching 
the first few pseudopotentials for the two problems. This 
form is motivated by the following observations: (i) The 
Fourier transform of the exact effective interaction for 

w ^ isi^i^ V^^{q) = — ( 1 — ^ ) • Its inverse Fourier 

transform V^^''{r) = ;^ + pr + 475-, however, is ill-behaved 
in that it yields divergent pseudopotentials for relative 
angular momenta to = 0, 1. (ii) Regularizing the interac- 
tion as in Eq. ([T]) removes the short distance divergences 
without significantly altering the long-distance behavior, 
(iii) Adding short-range Gaussian terms and fitting the 
first few pseudopotentials takes care of the short range 
part of the interaction without affecting the long-distance 



behavior. The constants Aj are arbitrary; we choose 
A3 — 1, Ac, = 10, ^7 = 100 to maximize the efficiency of 
our calculation. The interaction in Eq. ([1]) reproduces all 
second LL pseudopotentials almost exactly. The trans- 
verse thickness w is modeled through a square quantum 
well potential with the electronic wave function given by 
ip{z) — \/^cos(^^^ in the transverse dimension. Table 
H] gives Bi, Ci for various values of transverse thickness w, 
as well as the greatest relative error in pseudopotentials 
due to the approximations made in the form of Eq. ([1]) . 
The next step is to determine the interaction 
pseudopotentials^^ for composite fermions in the second 
A-level following Refs. Il6ll22l . by evaluating the energy of 
the state with two composite fermions in relative angu- 
lar momentum m state, the wave function for which can 
be constructed explicitly according to the standard CF 
theory^. (Another formalism22, has also been used for 
a treatment of the inter-CF interactions at the Hartrce- 
Fock level. However, that approach is designed for long 
distance physics and is not reliable for absolute energy 
comparisons of competing CF states.) Our calculation 
is based on the Monte Carlo method in the spherical 
geometryJ^^, in which electrons move on the surface of a 
sphere and a radial magnetic field is produced by a mag- 
netic monopole of strength Q at the center. Here 2Q(j>o 
is the magnetic flux through the surface of the sphere; 
(f>o — hc/e; and 2Q is an integer according to Dirac's 
quantization condition. The single particle states are 
monopole harmonics^ YQirm where I = Q + n is the 
angular momentum with n = 0, 1, . . . being the LL in- 
dex, TO = — /, — Z -|- 1, . . . , Z is the z-component of angular 
momentum. Composite fermion states are defined as^i^^ 



^^^ = ^LLL X{{U^VJ - V.Ujf^, 



(2) 



i<j 



where u = cos(0/2) exp(— i(/)/2), v = sin(6'/2) exp(i0/2), 
$ is a Slater determinant of YqimS, and T'lll is the 
lowest LL projection. (The planar geometry equivalents 
are obtained, apart from the Gaussian factor e~ 4 ^^ '^'1 , 
by the substitution {uiVj ~ ViUj) <=^ {zi — Zj), 
where Zi = xi — iyi denotes the coordinates of the «th 
particle on the plane.) For the pseudopotentials V^'^^ 
in the second A level, we consider states that contain 
two composite fermions in the second A level above a 
fully occupied lowest A level. In the spherical geome- 
try the pseudopotentials are size, or N dependent. For 
any specific N the pseudopotential Vm is the energy 
of two composite fermions at relative angular momen- 
tum TO. The interaction (Eq. JT])) is evaluated for such 
states assuming that r in Eq. 1 is the chord distance. 
To allow for a comparison between systems with differ- 
ent numbers of particles, an additive constant is cho- 
sen to fit the largest to pseudopotential to the expected 
asymptotic value Vm = '^~^ 2rrrn+T) between point-like 
charge e/3 objects at long distances. (The prefactor ac- 
counts for both the fractional charge and the conversion 



factor for energies expressed in terms of 



and 



££• 



wie.B 


B3 


Bs 


Br 


Co 


Ci 


Ci 


Ca 


rel. error at m - 


= 9 





1 


2.25 





-29.6652 


25.9333 


-5.7924 


0.35502 


-5 X 10"^ 




0.6 


0.98824 


2.1447 


-0.64764 


-26.7113 


22.8682 


-5.05705 


0.308495 


-4 X 10-*^ 




1 


0.96733 


1.96027 


-1.73059 


-23.7501 


20.1878 


-4.46304 


0.27412 


-3 X 10-"^ 




1.2 


0.95295 


1.83554 


-2.42498 


-22.4375 


19.2136 


-4.2825 


0.266517 


-1 X lO"*^ 




1.8 


0.89414 


1.3427 


-4.85138 


-19.8480 


18.4978 


-4.41921 


0.299473 


7 X 10-*^ 




2 


0.86931 


1.14304 


-5.68157 


-19.53315 


19.0897 


-4.72306 


0.332065 


1 X 10"^ 




2.4 


0.81181 


0.69985 


-7.17964 


-19.7498 


21.4679 


-5.70842 


0.429353 


2 X 10"* 




3 


0.70595 
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TABLE I: The coefficients of the second Landau level effective interaction (Eq. ([T])) as a function of quantum weU thickness. 



where ig = V^is is the magnetic length for composite 
fermionsi^ and is = \Jhc/eB.) Then V^^ is obtained 

from an linear extrapolation to the N —^ oo limit of Vm 
as a function of 1/iV. We used 18 < A^ < 50 for this ex- 
trapolation. As seen in Fig. [U the inter-CF interaction 
smoothly connects to the interaction between point-like 
charge e/3 objects at long distances, with significant de- 
viation at short distances. The assymptotic expression 
is used for m > 17, where the numerically calculated 
pseudopotential is fitted to the assymptotics modulo an 
additive constant. As expected, the transverse thickness 
weakens the short-range part of the interaction. 
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FIG. 1; (Color online) The pseudopotentials for two compos- 
ite fermion quasiparticles in the n — 1 LL, and the asymptotic 
form {2n + l)''^/'^V^^ = 3"^^^ 'arTm+i? ■ The solid line on the 
asymptotic form is a guide to the eye. 



ated, for which we use the forroi^ 



V(-\r*)=iY: 



{w) *8i+4 — r* 



(2n+ 1)^5/2 ^ g 



\i=0 



2 

(3) 

where r* is the distance measured in units of i*g . With 
six parameters (Table |lTl , all odd pseudopotentials Vm 
from ?n = 1 to ?n = 13 can be fitted exactly, and Eq. ([3]) 
already has the correct long-range behavior. 
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Armed with the CF pseudopotentials, we finally map 
the system into that of fermions at i'' — ly* — 1 = 1/2 
in the lowest Landau level. To the pseudopotentials of 
Fig. ID) an effective real-space interaction can be associ- 



FIG. 2: (Color online) Cohesive energy for the stripe (top) 
and bubble (bottom) phases of composite fermion quasipar- 
ticles in the second A-level of second Landau level. 
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TABLE II: The coefficients of the effective inter-CF interaction in the second A level of the second Lnadau level (Eq. (|3])) as a 
function of quantum well thickness. 



III. VARIATIONAL STATES 

We first consider charge-density-wave states of com- 
posite fermions, both stripes and bubble crystals; analo- 
gous states have proved to be relevant for half-filled elec- 
tronic LLs with high LL inde»2i. (We note that liquid 
crystalline phases of electrons, possibly with nematic or- 
der, have also been considered in high Landau levels^; 
we do not consider in this work analogous phases for com- 
posite fermions.) In the Hartee-Fock scheme the cohesive 
energy of these states is^i 



(2^)3 



2NL.,Ly 



Y, C/HF(g)A(-q)A(q), (4) 



q#0 



which is defined as the interaction energy measured from 
the uniform Hartree-Fock state 



_ t/(g = 0) 
-C/O — t: '-' ■ 



(5) 



This expression is based on the assumption that the 
CF-background and the background-background inter- 
action also have the same form as the CF-CF inter- 
action in Eq. ([3]); because the first two are identi- 
cal for all uniform states, their actual form is not 
relevant to the energy comparisons, and can be cho- 
sen according to convenience. The quantity A(q) = 
^Efce~'"'"^^^^'(4+afc-) in Eq. g]) is the orbit-center 
density, and we define C/hf(9) = U{q) - {£*B)'^U{qi*g), 
U(q) = T/(»)(q)e-59'(^B)', and k± = k^qy/^. For the 
stripe phase this reduces to 



T-istripc 
^coh 



2iy'{il 



E, , 2 sm ^^ — 



g#0 



AsQ 



where q = -4^, and for the bubble crystal 



Tibubblc 



i?cT''"= = ^E^HF(<z) 
I/' ^ — ' 



R 



Ji(gi?) 



(6) 



(7) 



with R = Ki,J^/lv' /2Tr, A = (\/3/2)Ag, and q = nei 

-y- 



me2 with ei = -^y and 63 - f^ic - ^^y. The 



V^\b Aft ysAb 

parameters As and Ah denote the period of the stripe 
and bubble phases, respectively. The results shown in 
Fig. [21 indicate interesting differences from electrons in 
higher LLs^^ and also from composite fermions in the 
lowest Landau leveU^. The stripe phase is favored in an 
intermediate range 0.12 "^ v' < 0.4; the bubble crystal 
is better for v' < 0.12; and the two are very competitive 
close to half-filling 0.4 "^v'k 0.5. The periods are 8.5 < 
Ab < 10 and As « 8 — 9, apart from small u' where the 
stripe phase is irrelevant. 

The sharp short-range repulsion between composite 
fermions (Fig. [T|) suggests the possibility of further vor- 
tex attachment to account for correlations between them. 
We consider the PfafRan wave function^. 



^Pf = 



]^(MiWj - WjUj)^Pf 



K'j 



1 



ViU, 



(8) 



which describes an incompressible p-wave paired state of 
composite fermions, and the compressible CF Fermi sea. 



^ 



CFFS 



t'lll n(^ 



^$ 



FS- 



(9) 



i<3 



For a comparison with the charge-density- wave states we 
calculate the cohesive energy of VP^^ and ^^^^^ by sub- 
tracting the energy of the uniform Hartree-Fock state 
(Eq. (O). The cohesive energies of ^^^ and v]/Cffs 
are shown in Fig. [31 The extrapolation is based on 
iV = 30 — 100 particle systems for the PfafRan, and 
N — 36, 49, 64 for the Fermi sea. In the thermodynamical 
limit N ^ 00 both the PfafRan and the CFFS states are 
energetically favored over the charge-density-wave states 
f Table HIT)) . The PfafRan wave function has higher en- 
ergy for the whole range of transverse thickness studied, 
making PfafRan-CF pairing unlikely as a mechanism for 
incompressibility at 2 -I- 3/8. 
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FIG. 3: (Color online) Thermodynamic extrapolation of the 
cohesive energy of the Pfafhan and the CF Fermi sea wave 
functions. The solid line extends to the right up to the small- 
est system included the linear fitting (A'^ = 30 and 36, respec- 
tively) . 
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-0.03050(2 
-0.03102(2 
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TABLE III: The cohesive energy of the Pfaffian and the CF 
Fermi sea wave functions. 



IV. EXACT DIAGONALIZATION 

We have confirmed the above conclusion by performing 
exact diagonalization for fcrniions at z/' — 1/2 interacting 
with the potential given in Eq. ([3]). Exact diagonalization 
at 2Q = 2N — 3 for TV = 8, 10 shows the ground state 
to be nonuniform, consistent with the variational study 
ruling out the Pfaffian wave function. However, we find 
that the quantum number of the ground state at 2Q = 



2N — 2 agrees with the prediction of the CF theory for 
A^ = 4 — 11 (see Table HV)) . The same quantum numbers 
result for the half-filled lowest Landau level, where the 
CF Fermi sea state is well establishe d^°'^^ . 
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TABLE IV: The ground state angular momentum for 2Q = 
2N — 2 on the sphere, along with its interpretation. The state 
with n filled A levels is denoted by <!>„. 



V. COMPOSITE FERMION 
DIAGONALIZATION 

While the noninteracting CF Fermi sea is compressible, 
the possibility that the residual inter-CF interaction may 
give rise to incomprcssibility cannot be a priori excluded. 
An investigation of this physics requires larger systems 
than can be addressed in exact diagonalization. We have 
studied this problem by a perturbative process called CF 
diagonalizatiorti^ii^. (Notice that we form ^CFs out of 
^CFs, and our starting point is a residual inter-^CF in- 
teraction V^^.) The wave functions for noninteracting 
composite fermions, for ground as well as excited states, 
can be constructed by analogy with the system of non- 
interacting electrons at an effective filling-^. These are 
given by 



*Q = -Plll I I \u,v 



n(^ 



luj 



ViUjY^Q* 



(10) 



(here, $q* is the wave function for the ground or ex- 
cited state at Q* = 0), and form bands separated by 
an effective CF cyclotron energy. At the n-th order CF 
diagonalization, we diagonalize the residual inter-CF in- 
teraction in the truncated space of correlated wave func- 
tions of the lowest n + \ bands, using the Metropolis 
Monte Carlo method^'^. It is necessary to go to at least 
second order to allow hybridization of the uniform (with 
orbital angular momentum L = 0) state with "excited" 
bands. At the second order, the ground state has L = 0, 
with a very small gap which depends only weakly on the 
layer thickness, as shown in Fig. SI The enhancement 
of the gap from iV = 16 to A^ = 25 suggests possible 
establishment of incompressibility due to the residual in- 
teraction between composite fermions. While our data 
in the N < 25 range, unfortunately, do not allow for an 
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FIG. 4: (Color online) The collective mode gap at ;/ = 2 + 
3/8 from second-order CF diagonalization for A^ = 9, 16, 25 
particles. The bars indicate the statistical uncertainty arising 
from Monte Carlo sampling. 



VI. CONCLUSION 

In conclusion, modeling the system at i^ = 2 + 3/8 
as filling factor 3/2 of fully spin polarized composite 
fermions in the second electronic Landau level, we have 
considered many possible structures by several meth- 
ods. Our study suggests the possibility of a very delicate 
FQHE here due to residual interactions between compos- 
ite fermions, but with a state distinct from the Pfaffian 
state. At the end, we note that even though our CF di- 
agonalization approach gives the low energy spectrum, it 
does not provide a simple wave function for the ground 
state, which has often been very useful in achieving a 
physical understanding of the physics of a FQHE state. 
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